Nonlinear dynamics and surface diffusion of diatomic molecules 
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The motion of molecules on solid surfaces is of interest for technological applications, but it is 
also a theoretical challenge. We study the deterministic and thermal diffusive dynamics of a dimer 
moving on a periodic substrate. The deterministic motion of the dimer displays strongly nonlinear 
features and chaotic behavior. The dimer thermal diffusive dynamics deviates from simple Arrhenius 
behavior, due to the coupling between vibrational and translational degrees of freedom. In the 
low-temperature limit the dimer diffusion can become orders of magnitude larger than that of a 
single atom, as also found experimentally. The relation between chaotic deterministic dynamics and 
stochastic thermal diffusion is discussed. 
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I. INTRODUCTION 

The surface diffusion of single adatoms has been in- 
tensively studied over the last decades, [1-3] due to its 
importance in thin film and crystal growth. [4] Once 
individual atoms are adsorbed on a surface they can 
meet, thus forming larger clusters. However, the diffu- 
sion of even the simplest cluster, a dimer, on a surface 
is by far not yet understood. [5-13] The diffusion dy- 
namics can be strongly affected by the coupling of the 
intramolecular motion to the translational motion of the 
centre of mass (CM) of the cluster. [12-15] Herein, we 
present a simple, one-dimensional model for studying the 
Hamiltonian and diffusive dynamics of a dimer, which 
is relevant to systems where quasi-one-dimensional mo- 
tion takes place. [16] The deterministic dynamics of this 
model is characterized by a complex behavior, dominated 
by non-linear effects, parametric resonances and chaotic 
features. At variance with the case of a single atom, at 
T ^ the role of the internal degrees of freedom of the 
dimer is responsible for deviations from activated behav- 
ior of the diffusion coefficient. In Section II we briefly 
outline our model. In Sections III and IV we discuss 
the nonlinear deterministic and thermal dynamics, re- 
spectively, and compare the two situations in Section V. 
Concluding remarks are given in Section VI. 

II. MODEL 

We consider the deterministic and thermal dynamics of 
a dimer moving on a periodic one-dimensional substrate. 
The particle-substrate interaction is a sinusoidal function 
of amplitude 2Uq and period a, and the interparticle in- 
teraction is given by a harmonic potential with spring 
constant K and equilibrium length I. We use Langevin 



dynamics to deal with finite temperature T. The equa- 
tions of motion for the two atoms of mass m and of co- 
ordinates x\ and X2 composing the dimer are given by 
Equation (1): 
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where the effect of finite temperature T is taken into 
account by the stochastically fluctuating forces fi, satis- 
fying the conditions < fi(t) >= and < fi(t)fj(0) >= 
2mrjkBT6ijS(t), and by the damping term mrjii. In the 
following, we will use representative values of the parame- 
ters, a = 0.25 nm, U = 0.2 eV, m = 5 x 10~ 26 kg, 77 = 0.7 
ps _1 . The values of K and I will be given in the caption 
of each figure. We performed molecular dynamics (MD) 
simulations, integrating the equations of motion using a 
velocity- Verlet algorithm, with a time step A = 10~ 16 s 
and averaging the trajectories over several thousands of 
realizations in the case of thermal diffusion, in order to 
reduce the statistical noise. 



III. HAMILTONIAN DYNAMICS 

First we consider the Hamiltonian dynamics (ft = 
and r\ = in Equation (1)). It is convenient to rewrite 
Equation (1) in terms of the CM coordinate xcm = 
(x\ + X2)/2 and of the deviations from equilibrium of 
the internal coordinate x r — X2 — x\ — I, obtaining Equa- 
tion (2): 



mxcM 
m'x r 



_ 2-itUq 



sin 



= -2Kx r 



«u 



cos 



cos I 

( 2-kxcm 



liXr+l)} 

' Sin [f ( Xr + I)] 

(2) 



1 



We have considered the case of a commensurate dimcr 
(I = a) starting at equilibrium with a given initial kinetic 
energy E^ in , a case which allows some analytical results 
for the initial phase of the motion, showing the role of in- 
ternal vibrations on the dynamics, to be obtained. [12] In 
fact, for a rigid dimer with xi(0) = ±2(0) = Vo, the mini- 
mum kinetic energy for the CM to overcome the p otentia l 
barrier is E® in = mv\ = AUq- Hence, for u < y/AUo/m, 
the motion of the CM is oscillatory, while a drift regime is 
attained for v > y/AU^/ra. Conversely, when the dimer 
is allowed to vibrate, the coupling between the CM and 
the internal motion makes it possible for the CM of the 
dimer to overcome the potential barrier AUq for values 
of vq below the threshold y/4Uo/m. In fact, if the in- 
ternal motion is excited, it can happen that one particle 
remains in the minimum and the other reaches the near- 
est maximum. From the energy balance, Equation (3) 
follows: 

El in = 2U + l -K{a/2f (3) 

If K is sufficiently small, the right-hand side of Equa- 
tion (3) can be smaller than 4£/ . This is the situation 
shown in Figure 1 , where the CM motion is rather irreg- 
ular, behaving in a chaotic fashion. 



IV. THERMAL DIFFUSION 

The thermal diffusive behavior of the dimer is char- 
acterized by computing the diffusion coefficient D from 
the mean square displacement < x^, M {t) >, as in Equa- 
tion (4): 

D = lim <X cm^ > (4) 

For temperatures that are small compared to the energy 
barrier, Kramer's theory [17] predicts an activated Ar- 
rhenius behavior of diffusion, given by Equation (5): 

D = D Q cxp(-E a /k B T), (5) 

where the activation energy E a and the prefactor D do 
not depend on T. In the limit of vanishing energy bar- 
rier, the adatom diffusion coefficient obeys the Einstein 
relation and is given by D = fcgT '/ '(mrf), twice the value 
of the dimer D = ksT '/ '(2mr)) . The diffusion coefficient 
of a single adatom shown in Figure 2 has an activation 
energy corresponding to the energy barrier 2Uq — 0.4 eV, 
except at very high temperatures owing to finite barrier 
effects. [18] 




FIG. 1. Dynamics given by Equation (2) for K = 0.2 Nm _1 
and vo = 380 ms _1 , fi — and 77 = 0. (a) The CM motion 
is plotted; (b) the deviations from equilibrium of the internal 
coordinate are plotted. 
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FIG. 2. Diffusion coefficient D as a function of l/(fcsT) 
for the adatom and the dimers with different values of / and 
K = 2 NnT 1 . The points are the result of the simulations 
and the lines represent fits to the data in the low-temperature 
regime. 



This chaotic regime occurs for weakly bound dimers in a 
velocity window around the threshold \fAUo, and can be 
characterized by Lyapunov exponents and power spec- 
tra. [12] We have also shown that, for larger values of 
K, the drift motion of the CM may excite the internal 
vibrations by a parametric resonance in a velocity win- 
dow around t wice the natural stretching frequency of the 
dimcr loq = y/2K/m. 



For a rigid dimer, E a — AUq for I = a down to E a = for 
I = a/2. If the dimer is not rigid, the activation energy 
is a non trivial function of the equilibrium length and 
clastic constant. In Figure 2, we show the diffusion coef- 
ficient for non rigid dimers with different values of /. We 
find E a = 0.6 eV for I = a, E a = 0.3 cV for I = a/2 and 
E a = 0.34 eV for the incommensurate case I — r g a, where 
T g = (1 + \/5)/2 is the golden mean. The finite-barrier 
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corrections to the Arrhcnius behavior at high temper- 
atures (ksT > 2Uo) are rather pronounced for / =/= a. 
Deviations from the Arrhcnius law due to incommensu- 
rability are also found for larger clusters. [19] 

It has been suggested that non- Arrhcnius behavior can 
result also from dynamical effects related to the internal 
motion of the dimcr. [13,15] The role of the internal vi- 
brations on the diffusive behavior is illustrated in Fig- 
ure 3(a), where we compare the diffusion coefficient of 
the rigid and non rigid incommensurate dimer (I = r g a). 
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FIG. 3. 



(a) Diffusion coefficient as a function of l/(fcsT) 
for the incommensurate dimer with I = r g a and K — 2 Nm -1 . 
The intramolecular length was kept fixed or not fixed, as in- 
dicated by the labels. Fitted activation energies [eV] are also 
reported, (b) Dynamical equilibrium length < / > t as a func- 
tion of l/(k B T). 

It is clear that we can define a unique value of E a for 
the rigid dimer, whereas the activation energy is, in gen- 
eral, temperature dependent when the dimer is allowed 
to vibrate. The temperature dependence of E a is linked 
to a temperature dependent misfit < I > t , induced by 
the dynamics, as shown in Figure 3(b). 

Moreover, the dimer can even diffuse faster than the 
adatom, at least for I ^ a and low temperatures. En- 
hanced diffusivity of dimers and small clusters is also 
found by theoretical studies of ID diffusion in molecular 
sieves [20] and in zeolite crystals. [21] Mitsui et al. [22] 
measured water diffusion on Pd(lll) at low temperature 
(T ~ 40 K), finding the mobility of dimers and larger 
clusters to be 3 — 4 orders of magnitude larger than that 
of adatoms. This experimental result is compatible with 
our findings. In fact, extrapolation to low temperature 
of our results of Figure 2 shows that the dimer diffusion 
for I 7^ a can be orders of magnitude higher than for the 
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V. RELATION BETWEEN DETERMINISTIC 
AND THERMAL DIFFUSION 

We find that the chaotic dynamics discussed in Sec- 
tion III can give rise to a diffusive behavior, even in the 
absence of thermal fluctuations. The role of the heat 
bath is played by the exchange between translational and 
internal motion which, owing to the nonlinearity of the 
system, can occur in a random manner. Figure 4 shows 
a comparison between the deterministic and the thermal 
mean square displacements. 
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FIG. 4. Mean-square displacement for (a) the deterministic 
and (b) the thermal motion for I = a (thick solid lines), (a) 
K = 0.1 Nm" 1 , v = 475 ms _1 , r) = 0. (b) K = 0.2 Nm" 1 . 
The thin solid lines are linear fits for large t; the dotted lines 
are power-law fits with exponent ~ 1.5 for small t. The data 
plotted in (b) were obtained by averaging the trajectory over 
3000 realizations. 

In the deterministic case < x 2 CM > represents a time av- 
eraging taken by displacing the time origin, [23] whereas 
at T ^ it is an average over realizations. The long-time 
behavior is linear (diffusive) even for the Hamiltonian dy- 
namics in the chaotic regime. In both cases we also find 
a transient superdiffusive regime (< x 2 CM >oc t a , with 
a ~ 1.5). Anomalous diffusion with a ~ 7/5 has also 
been observed in a model of adatom surface diffusion in 
two dimensions. [24] Anomalous diffusion might be re- 
lated to long jumps of the dimer, that is to trajectories 
that move over multiple surface minima, [24,25] which do 
occur during our simulations. Furthermore, it has been 
claimed that a deterministic diffusive behavior leads to a 
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non-Arrhenius dependence of the thermal diffusion coef- 
ficient, [26] which is compatible with our results. 



VI. CONCLUSIONS 

We have presented a one-dimensional model to de- 
scribe the diffusive dynamics of dimers on periodic sur- 
faces. We have shown that the coupling between trans- 
lational and vibrational degrees of freedom can lead to 
Hamiltonian chaotic motion and to non-Arrhenius be- 
havior of thermal diffusion. We have also pointed out 
the relation between deterministic and thermal diffusion, 
ft would be interesting to enhance the complexity of the 
system by considering orientational degrees of freedom 
in two dimensions and anharmonic intramolecular po- 
tentials. 



ACKNOWLEDGEMENTS 

This work was supported by the Stichting Funda- 
menteel Onderzoek der Materie (FOM) with financial 
support from the Nederlandse Organisatie voor Weten- 
schappclijk Onderzoek (NWO), project 0f5.000.031. The 
authors wish to thank Ted Janssen for contributing to 
this work and Sergey Krylov for useful and stimulating 
discussions. 



[1] H. Risken, The Fokker-Planck Equation, Springer, Berlin, 

1989, chapter ff. 
[2] G. L. Kellogg, Surf. Sci. Rep. 1994, 21, 1. 
[3] T. Ala-Nissila, R. Ferrando, S. C. Ying, Adv. Phys. 2002, 

51, 949. 



[4] J. A. Venables, G. D. T. Spiller, M. Hanbucken, Rep. 

Progr. Phys. 1984, 47, 399. 
[5] R. Wang, K. A. Fichthorn, Phys. Rev. B 1993, 48, 18288. 
[6] a) O. M. Braun, Surf. Sci. 1990, 23, 262; b) O. M. Braun, 

Phys. Rev. E 2001, 63, 011102. 
[7] a) F. Montalenti, R. Ferrando, Phys. Rev. Lett. 1999, 

82, 1498; b) F. Montalenti, R. Ferrando, Phys. Rev. E 

2000, 61, 3411. 

[8] G. Boisvert, L. J. Lewis, Phys. Rev. B 1997, 56, 7643. 
[9] A. S. Kovalev, A. I. Landau, Low Temp. Phys. 2002, 28, 
423. 

[10] A. H. Romero, A. M. Lacasta and J. M. Sancho, Phys. 
Rev. £2004, 69, 051105. 

[11] S. Gongalves, V. M. Kenkre and A. R. Bishop, Phys. Rev. 
B 2004, 70, 195415. 

[12] C. Fusco, A. Fasolino, Eur. Phys. J. B 2003, 31, 95. 

[13] C. Fusco, A. Fasolino, This Solid Films 2003, 428, 34. 

[14] P. Deltour, J.-L. Barrat, P. Jensen, Phys. Rev. Lett. 
1997, 78, 4597. 

[15] S. Yu. Krylov, Phys. Rev. Lett. 1999, 83, 4602. 

[16] a) U. Kiirpick, Phys. Rev. B 2001, 63, 045409; P. J. 
Feibelman, Phys. Rev. B 2000, 61, R2452; b) F. Mon- 
talenti, R. Ferrando, Surf. Sci. 1999, 432, 27. 

[17] P. Hanggi, P. Talkner, M. Borkovec, Rev. Mod. Phys. 
1990, 62, 251. 

[18] F. Montalenti, R. Ferrando, Phys. Rev. B 1999, 59, 5881. 
[19] J. C. Hamilton, Phys. Rev. Lett. 1996, 77, 885. 
[20] D. S. Sholl, K. A. Fichthorn, Phys. Rev. Lett. 1997, 79, 
3569. 

[21] a) D. S. Sholl and C. K. Lee, J. Chem. Phys. 2000, 112, 
817; b) D. S. Sholl, Ind. Eng. Chem. Res. 2000, 39, 3737. 

[22] T. Mitsui, M. K. Rose, E. Fomin, D. F. Ogletree, M. 
Salmcron, Science 2002, 297, 1850. 

[23] M. P. Allen, D. J. Tildesley, Computer Simulation of Liq- 
uids, Clarendon Press, New York, 1988, pp. 185-187. 

[24] R. Guantes, J. L. Vega, S. Miret-Artes, Phys. Rev. B 

2001, 64, 245415. 

[25] a) D. S. Sholl and R. T. Skodje, Physica D 1994, 71, 168; 
b) D. Sholl and K. A. Fichthorn, Phys. Rev. E 1997, 55, 
7753. 

[26] D. I. Kopelevich, H.-C. Chang, Phys. Rev. Lett. 1999, 

83, 1590. 



4 



